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Howard W. Sams & Co., Inc., shall have no liability or responsibility to the pur- 
chaser or any other person or entity with respect to any liability, loss or damage 
caused or alleged to be caused directly or indirectly by this product, including 
but not limited to any interruption of service, loss of business and anticipatory 
profits or consequential damages resulting from the use or operation of this 
product. 

RETURN PRIVILEGE 

This product may be exchanged within ninety (90) days from date of original pur- 
chase if defective as related to manufacturing, labeling or packaging. Except 
for such replacement, the sale or subsequent use of this product is without war- 
ranty or liability. Upon return, the defect must be explained in writing, detailing 
the failure of the product. Please return it, shipping charges prepaid, in its origi- 
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4300 West 62nd Street 
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return postage. 



Printed in the United States of Aerica 



TRS-80 I1STR0CT101S 



These programs load with the standard TRS-80 
"CLOAD 1 command. 

They may be saved to disk if you have a disk 
system by using the standard TRS-80 
procedures outlined in your DOS manual. 

If you have any problem loading these 
programs from the cassette. We suggest that 
you re-check the volume setting of your 
recorder and adjust if necessary. If 
problems persist, turn the cassette over and 
load the programs from the other side. All 
programs are recorded on both sides of the 
cassette. 

These programs have been left unprotected so 
that you may make a backup copy or alter 
these programs to your specific needs. 
These programs are meant for YOUR use only. 
We suggest, that you make a backup copy at 
this time to protect your software 
investment. 
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PLOTTING GRAPHS FOR VIDEO DISPLAY 

This TAPE illustrates how it is possible to graph functions or dis- 
crete data points on your video display in the form of histograms, 
Cartesian, semilogarithmic, log-log, or polar plots. Most of these 
programs are primarily written as subroutines so that they may be 
easily used in conjunction with your own special programs. 

HISTOGRAMS 

One of the simplest types of graphs is the histogram, or bar 
graph, which is usually used to display frequency distributions. 
Required input variables are: the number of data points N, the 
horizontal and vertical values X(l) and Y(l), respectively, for each 
point, and the graph title T$. Fig. 1 illustrates a simple driver 
program that can be used with the PHISTGM subroutine. The 
program assumes that each X-axis value has only a single corre- 

100 'MAIN PROGRAM FOR PHISTGM 

101 CLS:DIM X(101), Y(101) 

102 INPUT "NUMBER OF DATA POINTS TO BE PLOTED ";N 

103 FOR 1=1 TO N 

104 PRINT "POINT";l; 

105 INPUT "X, Y ";X(I),Y(I) 

106 NEXT I 

107 INPUT "TITLE FOR GRAPH ";T$ 

108 GOSUB 5000 'PHISTGM SUBROUTINE 

109 END 

Fig. 1. Simple main program for PHISTGM. 

sponding Y-axis value. Consequently, only a maximum of 101 
individual bars may be plotted. The DIM statement in line 101 must 
be included, as the coordinate values are stored before they are 
plotted on the display. 

Example 1 

Plot a histogram of the result of measuring the values of a sample 
of 55, 10% resistors, all of which are supposed to have the same 
nominal value, 100 ohms. As a result of our sampling, we find that 



the values vary, with several resistors having like values, as listed 
in Table 1. 

Table 1. Data for Example 1 



Point No. 


Value 


Frequency 


1 


90 


1 


2 


92 


1 


3 


93 


2 


4 


94 


4 


5 


95 


3 


6 


96 


4 


7 


98 


7 


8 


99 


8 


g 


100 


11 


10 


101 


4 


11 


103 


2 


12 


104 


3 


13 


106 


3 


14 


110 


2 

Total = 55 



Using the above data, Fig. 2A shows the results of executing the 
main program of Fig. 1 and the PHISTGM subroutine as one pro- 



NUMBER OF DATA POINTS TO BE PLOTTED ? 14 



POINT 1 X, 


Y ? 


90,1 


POINT 2 X, 


Y ? 


92,1 


POINT 3 X, 


Y ? 


93,2 


POINT 4 X, 


Y ? 


94,4 


POINT 5 X, 


Y ? 


95,3 


POINT 6 X, 


Y ? 


96,4 


POINT 7 X, 


Y ? 


98,7 


POINT 8 X, 


Y ? 


99,8 


POINT 9 X, 


Y ? 


100,11 


POINT 10 X 


, Y 


' 101,4 


POINT 11 X 


, Y ' 


' 103,2 


POINT 12 X 


, Y 


? 104,3 


POINT 13 X 


, Y 


? 106,3 


POINT 14 X 


, Y 


? 110,2 


TITLE FOR 


GRAPH ^ 100 OHM 



RESISTOR SAMPLE 



MIN - MAX X VALUES ARE: 90 
MIN - MAX Y VALUES ARE: 1 
MIN.MAX X-AXIS SCALE: ? 80,120 
MIN.MAX Y-AXIS SCALE: ? 0,15 

Fig. 2A. Plot of histogram of resistors, (output results) 



110 

11 



gram (refer to the appendix to see how a program stored on tape 
can be merged with another program already in memory). After all 
the 14 individual points and the graph title are entered, the 
display informs us as to the current ranges of the two axes. The 
program then asks us how we want the axes scaled, such as hav- 
ing the horizontal axis range from 80 to 120 while the vertical 
axis varies from to 15. Because of the method in which the axes 




Fig, 2B. Plot of histogram of resistors, (video plot) 

are scaled, it is best to have the vertical axis (i.e., the "frequency" 
axis) incremented by a number that is easily divisible by five, 
since there are five major divisions, or tic marks, on this axis. 
When the ENTER key is pressed, the resultant histogram is drawn 
on the video display, as shown in Fig. 2B. 

The histogram then remains on the screen until the ENTER key is 
pressed, after which we are asked whether or not we are satisfied 
with how the axes were scaled, as shown by the new entries in 
Fig. 3A and the resultant graph in Fig. 3B. The process of changing 
the axis scales can be repeated as many times as desired. 

ANY CHANGES IN SCALE FACTORS (YES/NO) ? YES 
MIN - MAX X VALUES ARE; 80 120 

MIN - MAX Y VALUES ARE: 15 

MIN, MAX X-AXIS SCALE: ? 50,150 
MIN, MAX Y-AXIS SCALE: ? 0,15 

Fig. 3A. Histogram entries for new scale. 




Fig. 3B. Histogram of resistors with new scales, (video plot) 



CARTESIAN PLOTS 

The program for generating a Cartesian, or rectangular X-Y co- 
ordinate plot is similar to the histogram. The PCARTXY subroutine 
can be used with the plot driver program of Fig. 1, except that the 
GOSUB 5000 statement in line 108 is now replaced with GOSUB 
5050. 

Example 2 

Graph the following 6 discrete (X,Y) points: (0,0.032), (0.02,0.135), 
(0.04,0.187), (0.06,0.268), (0.08,0.359), and (0.1,0.435). Fig. 4A shows 
the results after this data is entered. The resultant Cartesian plot 
is drawn on video display, as shown in Fig. 4B. The X-axis varies 

NUMBER OF DATA POINTS TO BE PLOTTED ? 6 



POINT 1 X 


Y ? 0,0.032 


POINT 2 X 


Y ? 0.02,0.135 


POINT 3 X 


Y ? 0.04,0.187 


POINT 4 X 


Y ? 0.06,0.268 


POINT 5 X 


Y ? 0.08,0.359 


POINT 6 X 


Y ? 0.1,0.435 



TITLE FOR GRAPH ? EXAMPLE 2-2 DATA 

CURRENT MIN - MAX X VALUES ARE: .1 

CURRENT MIN - MAX Y VALUES ARE: .032 

DESIRED MIN, MAX X-AXIS SCALE: ? 0,.1 

DESIRED MIN.MAX. Y-AXIS SCALE: ? 0..45 

Fig. 4A. Resylts afjer data is entered, (output resylts) 



.435 




Fig, 4B. Number of data points plotted, (video plot) 

from to 0.1 as Y ranges from 0.032 to 0.435, based on our data 
points. The X-axis was then scaled from to 0.1 while the Y-axis 
was conveniently set to range from to 0.45. 

On the other hand, suppose we are interested in graphing con- 
tinuous mathematical functions, instead of discrete data points. 
The following example illustrates how this is accomplished. 

Example 3 

The transient response of a given network has the known time 
varying solution: 



y(t) 



2.3te- 



Fig. 5 lists the main program used to plot this function, which is 
defined in line 106 in terms of the variables y, and time, t (the inde- 
pendent variable). Line 107 transforms these two variables to the 
variable arrays Y(l) and X(l), respectively, which are required by 
the PCARTXY subroutine. Since most functions are single-valued, 
there will be 101 continuous points, which are recognized by the 
DIM statement in line 101. When this main program is run, we are 
first asked to select the minimum and maximum values of time 
over which this function is to be plotted, after which we are asked 
for a suitable title for our graph. 

For this example, we are interested in how this function appears 
as time varies from to 2,0 seconds. When executed, we find that 



100 'SAMPLE PROGRAM FOR PLOTTING A FUNCTION USING PCARTXY 

101 CLS:DIM X(101),Y(101) 

102 N=101:l=1 

103 INPUT"MINIMUM AND MAXIMUM VALUES OF T TO BE PLOTTED 

";A,B 

104 C=(B—A)/100 

105 FOR T=A TO B STEP C 

106 Y=2,3*T*EXP(— 5*T) 'SAMPLE FUNCTION 

107 X(I)=T:Y(I)=Y:I=I+1 

108 NEXT T 

109 INPUT'TITLE FOR GRAPH ";T$ 

110 GOSUB 5050 

111 END 

Fig. 5. Main program for plotting a function of time. 



y ranges from to .17 as t varies from to 2.0, as shown in Fig. 6A, 
while the resulting plot is shown in Fig. 6B. The graph remains on 
the screen until the ENTEFl key is pressed, after which we are 
asked whether or not we wish to change the scale of either one 
or both axes. 

MINIMUM AND MAXIMUM VALUES OF T TO BE PLOTTED ? 0,2 

TITLE FOR GRAPH ? Y = 2,3*T*EXP(-5*T) 

CURRENT MIN - MAX X VALUES ARE: 2 

CURRENT MIN - MAX Y VALUES ARE: .169225 

DESIRED MIN,MAX X-AXIS SCALE: ? 0,2 

DESIRED MIN.MAX Y-AXIS SCALE: ? 0,.2 

(A) Output results. 




(B) Video plot. 

Fig. 6. Minimum and maximum values of T plotted. 



Periodic, or circular functions, such as-those which have sine or 
cosine terms are also easily graphed on the display, as illustrated 
by the following example. 



100 'SAMPLE PROGRAM FOR PLOTTING A FUNCTION USING PCARTXY 

101 CLS:DlM X(101),Y(101),H(6) * 

102 N=101:l=1 

103 INPUT'ENTER MINIMUM AND MAXIMUM VALUES OF X TO BE 
PLOTTED ";A,B 

104 C=(B-A)/100 

105 FOR X=A TO B STEP C 

106 Y=1.273*SIN(X)4-.424*SIN(3*X) + .255*SIN(5*X) 'FUNCTION TO 
BE PLOTTED 

107 X(I)=X:Y(I)=Y:I= 

108 NEXT X 

109 INPUT'TITLE FOR GRAPH ";T$ 

110 GOSUB 5050 

111 END 

(A) Main program. 



-.424*SIN(3*X) + .255*SIN(5*X) 
=1+1 



MINIMUM AND MAXIMUM VALUES OF X TO BE PLOTTED ? 0,6.28 
TITLE FOR GRAPH ? FOURIER SERIES - SQUARE WAVE 



CURRENT MIN - MAX X VALUES ARE: 
CURRENT MIN - MAX Y VALUES ARE: -1.18677 
DESIRED MIN.MAX X-AXIS SCALE ? 0,6.28 
DESIRED MIN, MAX Y-AXIS SCALE ? -1.5,1.5 

(B) Output results. 



6.28 



1.18653 






« ar as as 



1.256 2.512 3.768 5.024 6.28 

FOURIER SERIES-SQUARE WAVE 
HIT (ENTER) TO CONTINUE?, 



(C) Video plot. 
Fig. 7. Plotting a function using PCARTXY, 



Example 4 

Plot the first three terms for the Fourier series of a square wave, 
such that: 

f(x) = 1.273sin(x) + 0.424sin(3x) + 0.255sin(5x) 

Fig. 7A lists the main program used, plotting the function in terms 
of radians. Fig. 7B shows the output results along with the plot of 
this function in Fig. 7C. 

SEMILOGARITHMIC PLOTS 

The PSEMILOG subroutine permits the plotting of points on a 
semilogarithmic coordinate system, which is logarithmic in X 
and linear in Y. Depending on the range of values for the hori- 
zontal axis, the subroutine automatically scales the X-axis for 
values with a minimum of 4 cycles. However, the user is able to 
set the range for the Y-axis scale. 

LOG-LOG PLOTS 

The LOGLOG subroutine is used for plotting points whose range 
requires logarithmic scales for both axes. The main program of 
Fig. 8 can be used to run with the subroutine LOGLOG. 

100 'MAIN PROGRAM FOR LOGLOG SUBROUTINE 

101 CLS:DIM X(101),Y(101),H(6),V(6) 

102 INPUT"NUMBER OF DATA POINTS TO BE PLOTTED ";N 

103 FOR 1=1 TO N 

104 PRINT"POINT";l; 

105 INPUT"X, Y ";X(I),Y(I) 

106 NEXT I 

107 GOSUB 5200 'LOGLOG SUBROUTINE 

108 END 

Fig. 8. Simple main program for LOGLOG. 

POLAR PLOTS 

Polar plots are frequently used in field theory and control systems, 
where the plot describes both a radial distance as well as the 
angle (as measured counterclockwise from the horizontal) asso- 
ciated with it. Examples are the electric field distribution of an 
antenna, or Nyquist plots of the stability of feedback systems. 



The VPOLAR program permits the plotting of either discrete data 
points or functions expressed in polar coordinates. The polar 
function must be defined in line 152. The variable I is in degrees 
and is converted to radians by dividing it by 57.29577. As exam- 
ples, Table 2 lists the polar forms of a yariety of frequently used 
functions, where the angle is expressed- in radians. 

Table 2. Common Polar Functions 



Name 


Equation 


Circle: 


r = A cos9 or r = |A sin9| 


Cardioid: 


r = A(l + cosS) or r = |A(1 + sin9)| 


Lemninscate: 


r = |ACos(n9)| or r = |Asin(n9)| 


Limacon: 


r = !A + Bcos9| or r = [A + Bsin8| 


N-Leaved Rose: 


r = t Acos(n8)| or r = Asin(n8)| 


Archimedian Spiral: 


r = A9 


Hyperbolic Spiral: 


r = A/8 


Logarithmic Spiral: 


r = e B8 



After the data points are entered or the polar function is deter- 
mined, the program asks for a graph title, and then determines the 
maximum radius that exists. We are then asked to enter a value 
for the radius which is the maximum scale for both axes. The 
plot is then shown on the video display, which contains markers 
30 degrees apart on the outer perimeter. In addition, the title, 
maximum radius, X-axis (DX) and Y-axis (DY) increments are also 
shown. 



Example 5 

Plot the cardioid 



r = 2(1 + cos(0)) 



so that there will be a point plotted every 4 degrees. Line 152 
should then read: 

152 R(K) = 2*(1+COS(l/57.29577)) 

When run, the maximum radius was found to be 4.0, giving the re- 
sults shown in Fig. 9. Other polar functions plotted using the 
VPOLAR routine are also shown in Fig. 9. 
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Fig. 9. Common polar functions using VPOLAR. (video plots) 
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APPENDIX 



MERGING A PROGRAM WlfH A SUBROUTINE 

A number of programs in this book were purposely written as sub- 
routines so that they may be made part of your own application 
programs. These include a variety of plotting routines, and mathe- 
matical functions, such as a four-quadrant arctangent function. 
Rather than retype the subroutine each time it is needed, one can 
easily merge, or combine, the subroutine which is already stored 
on tape with the main program which is presently in memory. 
To merge a subroutine already stored on tape with a main program 
in memory, perform the following steps: 

1. The line numbers of the subroutine must be larger than the 
highest line number of the main program. The subroutines 
in this book start with line numbers 5000 and higher, and no 
two subroutines have the same inclusive line numbers. 

2. While in BASIC ll's COMMAND mode, determine the values 
stored in locations 16633 and 16634 by typing: 

PRINT PEEK(16633),PEEK(16634) 

and then pressing the ENTER key. 

3. If the value stored in location 16633 is equal to or greater 
than 2, then execute the following two statements in the 
COMMAND mode: 

POKE 16548,PEEK(16633)-2 
POKE 16549,PEEK(16634) 

and then use the CLOAD command to load the program 
stored on cassette tape from the tape recorder, after which 
the following two statements are executed: 

POKE 16548,233 
POKE 16549,66 

The merged program can then be saved on tape using the 
CSAVE command. 
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4. If the value stored in location 16633 is less than 2, then exe- 
cute the following two statements: 

POKE 16548,PEEK(16633)+254 
POKE 16549,PEEK(16634)-1 

and then use the CLOAD command to load the program 
stored on cassette tape from the tape recorder, after which 
the following two statements are executed: 

POKE 16548,233 
POKE 16549,66 

The merged program can then be saved on tape using the 
CSAVE command. 
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SAMS SOS™ TECHNICAL SOFTWARE PROGRAMS 
FOR TRS-80, APPLE, AND OSI 

PLOTTING GRAPHS FOR LINE PRINTER (3 Programs) 

Cartesian Plots. Semi-Logarithmic Plots. Polar Plots. 

No. 26000-TRS-80 (tape); No. 26023-APPLE (tape); No. 26016-OSI (tape); 

No. 26009-APPLE (disk); No. 26030-OSI (disk). 

ACTIVE FILTER DESIGN (6 Programs) 

Low and High-Pass (Bessel, Butterworth, 1, 2, and 3-db Chebyshev). 

State-Variable Filter, Bandpass Filters with Q's less than 10 and 50, 

Staggered-Tuned Butterworth Bandpass Filters (2, 3, or 5 states). Notch Filter. 

No. 26001 -TRS-80 (tape); No. 26024-APPLE (tape); No. 26017-OSI (tape); 

No. 26010-APPLE (disk); No. 26031 -OSI (disk), 

DESCRIPTIVE STATISTICS & REGRESSION ANALYSIS (3 Programs) 

Descriptive Statistics (mean, standard deviation, variance, kurtosis, z-scores). 
Curvilinear Regression (linear, inverse, polynomial, exponential, logarithmic), 

Multivariable Linear Regression, 

No. 26002-TRS-80 (tape); No. 26025-APPLE (tape); No. 26018-OSI (tape); 

No. 26011 -APPLE (disk); No. 26032-OSI (disk). 

ELECTRONICS I (5 Programs) 

Zener Diode Voltage Regulator Design, 

555 Timer Design (monostable and astable circuits). Transistor Bias Parameters, 

SIngle-Stage Transistor Amplifier Design, Heat Sink Selection and Design. 

No. 26003-TRS-80 (tape); No. 26026-APPLE (tape); No. 26019-OSI (tape); 

No. 26012-APPLE (disk); No. 26033-OSI (disk). 

ELECTRONICS II (7 Programs) 

4-Quadrant Arctangent Function, 

Rectangular/Polar Conversion and Complex Number Mathematics. 

Minimum and Maximum Values of an Array. 

Roots of Polynomials with Real Coefficients. 

Inverse Laplace Transforms of a Transfer Function. 

Solution of Simultaneous Equations with Real and Complex Coefficients. 

No. 26004-TRS-80 (tape); No. 26027-APPLE (tape); No. 26020-OSI (tape); 

No. 26013-APPLE (disk); No. 26034-OSI (disk). 

ELECTRONICS III (8 Programs) 

Average and RMS Values of a Periodic Function. 

Fourier Series Expansion of a Periodic Function, 

Fourier Transform and Spectrum Plot, Analysis of Damped Oscillations. 

Impedance Matching Pads, PI-TEE (delta-wye) Transforms. 

No. 26005-TRS-80 (tape); No. 26028-APPLE (tape); No. 26021 -OSI (tape); 

No. 26014-APPLE (disk); No. 26035-OSI (disk). 

PLOTTING GRAPHS FOR VIDEO DISPLAY (5 Programs) 

Histograms. Cartesian Plots. Semi-Logarithmic Plots. Log-Log Plots. Polar Plots. 

No. 26006-TRS-80 (tape); No. 26029-APPLE (tape); No. 26022-OSI (tape); 

No. 26015-APPLE (disk); No. 26036-OSI (disk). 

Howard UU. Sams & Co., Inc. 

4300 W. 62nd Street, Indianapolis, Indiana 46268 USA 



